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Abstract. 

A generalized lens equation for weak gravitational fields in Schwarzschild metric 
and valid for finite distances of source and observer from the light deflecting body 
is suggested. The magnitude of neglected terms in the generalized lens equation is 

15 7T TTb^ 

estimated to be smaller than or equal to rr, where m is the Schwarzschild radius 

4 d' 2 

of massive body and d ' is Chandrasekhar's impact parameter. The main applications of 
this generalized lens equation are extreme astrometrical configurations, where Standard 
post-Newtonian approach as well as Classical lens equation cannot be applied. It is 
shown that in the appropriate limits the proposed lens equation yields the known 
post-Newtonian terms, 'enhanced' post-post-Newtonian terms and the Classical lens 
equation, thus provides a link between these both essential approaches for determining 
the light deflection. 



PACS numbers: 95.10.Jk, 95.10. Ce, 95.30. Sf, 04.25.Nx, 04.80.Cc 
1. Introduction 

Todays astrometry necessitates theoretical predictions of light-deflection by massive 
bodies on microarcsecond (^xas) level, e.g. astrometric missions SIM (NASA) or GAIA 
(ESA). In principle, an astrometric precision on microarcsecond level can be achieved 
by numerical integration of geodesic equation of light-propagation. On the other side, 
modern astrometric missions like GAIA determine the positions and proper motions of 
approximately one billion objects, each of which is observed about one hundred times. 
The data reduction of such huge amount of observations implies the need of analytical 
solutions, because the numerical investigation of geodesic equation is by far too time- 
consuming. 

The metric of a massive body can be expanded in terms of multipoles, i.e. monopole 
term, quadrupole term and higher multipoles [T71 Usually, the largest contributions 
of light-deflection originates from the spherically symmetric part (Schwarzschild) of the 
massive body under consideration. The exact analytical solution of light-propagation 
in Schwarzschild metric [7] inherits elliptic integrals, but their evaluation becomes 
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comparable with the time effort needed for a numerical integration of geodesic equation. 
Thus, approximative analytical solutions valid on microarcsecond level of accuracy are 
indispensible for a highly time-efficient data reduction. 

In the same way, exact lens equations of light-deflection have been obtained in 
[HI HH E]. However, such exact relations are also given in terms of elliptic integrals. 
Therefore, approximations of these exact solutions are also needed for a time-efficient 
data reduction. An excellent overview of such approximative lens equations has recently 
been presented in [3]. 

Basically, two essential approximative approaches are known in order to determine 
the light-deflection in weak gravitational fields: 

The first one is the standard parameterized post-Newtonian approach (PPN) [131 E] 
which is of the order O (m) . During the last decades, it has been the common 
understanding that the higher order terms O (m 2 ) are negligible even on microarcsecond 
level, except for observations in the vicinity of the Sun. However, recent investigations 
[T5l (H [HJ [16] have revealed that the post-post-Newtonian approximation [5] [6] , which 
is of the order 0(m 2 ), is needed for such high accuracy. Both approximations are 
applicable for d 3> m, where d being the impact parameter of the unperturbed light ray. 

The second one is the standard weak-field approximative lens equation, which is 
usually called the classical lens equation, see Eq. (67) in [9] of Eq. (24) in [I]. One 
decisive advantage of classical lens equation is its validity for arbitrarily small values of 
impact parameter d. The classical lens equation is valid for astrometrical configurations 
where source and observer are far ernough from the lens, especially in case of a d 
and b ^> d, where a = k ■ X\ and b = — fc • x , where x and X\ are the three- vectors 
from the center of the massive body to the source and observer, respectively, and k is 
the unit vector from the source to the observer. However, the classical lens equation is 
not applicable for determine the light-deflection of moons of giant planets in the solar 
system, because astrometrical configurations with 6 = are possible. 

Moreover, there are astrometric configurations where neither the standard post- 
Newtonian approach nor the classical lens equation are applicable, for instance binary 
systems. In order to investigate the light-deflection in such systems a link between these 
both approaches is needed. Such a link can be provided by a generalized lens equation 
which, in the appropriate limits, coincides with standard post-Newtonian approach 
and classical lens equation. Accordingly, the aim of our investigation is an analytical 
expression for the generalized lens equation having a form very similar to the classical 
lens equation. We formulate the following conditions under which our generalized lens 
equation should be applicable: 

1. valid for d = , a = X\ ^> m , b = x 3> m, 

2. valid for a = , d^> m , b ^ 0, 

3. valid for b = , d^> m , 

These conditions imply that the light-path is always far enough from the lens, thus 
inherit weak gravitational fields, i.e. small light deflection angles. In order to control 
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the numerical accuracy, the generalized lens equation is compared with the numerical 
solution of exact geodesic equation in the Schwarzschild metric (throughout the paper, 
we work in harmonic gauge): 

1 ~~ a - n 

900 — — -7—. > 9i0 — u ; 

1 + a 

gij = (1 + a) 2 Sij + ^ + a x l x J . (1) 



Here, a = — and m = — 5— is the Schwarzschild radius and M is the mass of the 
x & 

light-deflecting body, G is Newtonian constant of gravitation and c is the speed of light. 
Latin indices take values 1,2,3, and the Euclidean metric 5^ = 1(0) for i = j (i 7^ j). 



The absolute value of a three- vector is denoted by x = \x\ = J x\ + x\ + x\. The exact 
geodesic equation in Schwarzschild metric reads, cf. [H] 

2 — a fx ■ x 



x = — 



—c 2 - r= — x ■ x + a- 



„ a 1 — a , ... 

x + 2 ^i ~{x-x)x, 2 

x z 1 — ar 



1 + a) 3 1 — a 2 V x 

where a dot denotes time derivative in respect to the coordinate time t, and x is 
the three-vector pointing from the center of mass of the massive body to the photon 
trajectory at time moment t. The scalar product of two three-vectors with respect to 
Euclidean metric 5^ is a ■ b — dij a l V . The numerical solution of this equation will be 
used in order to determine the accuracy of approximative solutions. We abbreviate the 
angle between two three- vectors a and b by S(a, b), which can be computed by means 

of 8(a, b) = arccos — — . 

ab 

The paper is organized as follows: In Section |2] the post-Newtonian approach is 
presented. In Section [3] the steps of post-post-Newtonian approach relevant for this 
investigation are shown, and I will briefly summarize the main results of our article [11] . 
The generalized lens equation is obtained in Section H] and discussed in Section El A 
summary is given in Section |6j 



2. Post-Newtonian approximation 

Let us consider the trajectory of a light-signal in post-Newtonian Schwarzschild metric: 
g 00 = - 1 + 2 a + O (V 4 ) , g i0 = , 

9ij =5 ij + 2 1 a5 ij + 0(c- A ) . (3) 

Here, 7 is the parameter of the Parametrized Post-Newtonian (PPN) formalism, which 
characterizes possible deviation of the physical reality from general relativity theory 
where 7 = 1. The light-ray is being emitted at a position Xq at time moment t and 
received at position x\ at a time moment t±, see Figured! 

Light propagation is governed by the geodesic equation, in post-Newtonian order given 
by 

*= -(l + l)c 2 ^ + 2(l +1 ) a±i± 2 X) +0(c- 2 ). (4) 
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Figure 1. A geometrical representation of the boundary problem under consideration 
for a light-propagation from the source to the observer. 



The unit tangent vector at the point of observation is n = ' 7-^- , and the unit tangent 

|ac(ti)| 

R 

vector k — — , where R = x\ — xq and the absolute value is R = \R\. Furthermore, we 
R 

x(t) 

define the unit tangent vector at remote past: cr = lim . 

t— > — OO Q 

Up to post-Newtonian order, the differential equation can be solved analytically. 

The solution for the transformation between n and k reads 

u (1 i d' x x 1 -x Q -x 1 / 2 x 
„ = fc-(l + T )__ — + 0(m), (5) 

in terms of the coordinate- independent impact vector d', cf. Eq. (57) of 

d' = lim cr x (x(t) x cr) . (6) 

t— ^— 00 

This impact parameter is identical to Chandrasekhar's impact parameter [TT| [12]. that 

means in vectorial form d' = where -L is the orbital three- momentum and E is the 

energy of the photon on the light-trajectory; cf. Eq. (215) in chapter 20 of [7j. 

By means of sin</? = \n x fe|, we find the light-deflection angle tp = 5(n,k) in 

post-Newtonian approximation: 

m x xi - x ■ x x / 2 x 
V,= (l + 7)^ ^ + 0(m). (7) 

at ! xqX\ — xq-x\ Am 

Note that < 2, and therefore 99 < — — -. One problem of post-Newtonian 

Rx\ a' 

solution (J5J) or ^ is, that one can only state that the neglected terms are of order 
O (m 2 ), but their magnitude remains unclear. In order to make a statement about the 
upper magnitude of the higher order terms one needs to consider the geodesic equation 
in post-post-Newtonian approximation. 

3. Post-post Newtonian approximation 

Now we will consider the trajectory of a light-signal in post-post-Newtonian 
Schwarzschild metric: 

2oo= -l + 2a-2f3a 2 + 0(c- 6 ) , g i0 = , 
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gij = Stj + 2 7 a 5ij + e Uj + ^ j a 2 + O (V 6 ) . (8) 

The geodesic equation of light-propagation in post-post-Newtonian approximation is 
given by [TTJ 



, ax . ax (x • x) 

'— + 2 (! + 7) \ 



sc = - (1 + 7 ) c 2 — + 2(1+7) 



+ 2c 2 (/3- e + 2 7 (l + 7)) — + 2e 



+ 2(2(1 -/3) + e- 2 7 2 ) a2±( t" a:) +0(c- 4 ). (9) 



x 2 

The parameters /?, 7 and e characterize possible deviation of physical reality from general 
relativity theory (in general relativity /3 = 7 = e = 1). The solution of (j9]) and the 
transformation between the unit vectors n and k in post-post-Newtonian order has 
been given in [11], cf. Eqs. (108) and (109) ibid., and reads 

n = fc - (1 + 7) ^ + W)- (10) 

^72 j can be estimated to be smaller than or equal to 
- — 2. From Eq. (Tl0|) we obtain the expression 

^ =(1+7) ^ — ^ — + W 2 )- { ] 

The solutions ({TO]) and (fTT]) are identical to the post-Newtonian solution (jSJ) and (J?)), 
respecticvely. This fact means that the post-post-Newtonian terms in the metric (jSD and 
also the post-post-Newtonian terms in the geodesic equation (Q contribute only terms 

15 7T TTt^ 

which can be estimated to be smaller than or equal to - Therefore, the only 

difference between fTTTj) and (JTj) here is, that the post-post-Newtonian approximation 
allows to make a statement about the upper magnitude of the regular post-post- 
Newtonian terms. 



4. Generalized lens equation 

Usually, in practical astrometry the position of observer x\ and the position of light- 
deflecting body is known (here, the center of massive body coincides with the coordinate 
center), but the impact parameter d' is not accessible. Therefore, the solutions ( flDl or 
(TTTj) are not applicable in the form presented. Instead, one has to rewrite these solutions 
in terms of the impact vector of the unperturbed light ray 

d = kx( Xl xk). (12) 

For that one needs a relation between impact vector d' defined in Eq. ([6]) and impact 
vector d defined in Eq. ( fl2l) . Such a relation has been given in [TT] . cf. (62) ibid., and 
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reads (note, d ' — d + O (m) ) : 

d = d + (1 + 7) — — + O (m J . (13) 

Eq. (TT31) represents an quadratic equation for d', and these both solutions correspond 
to the two possible light-trajectories. A comparison of (TT3"]) with (TTTT) yields the relation 

d , = d + x 1 <p+ X0 + X ^~ R x 1 <p + O(m 2 ) , (14) 

where (p is given by Eq. (TTTT) and we have separated a term X ° 1 — — x\ ip which can 

R 

(m 2 \ 
-j-p \ • By inserting 

(114p into ( TTTT) we obtain an quadratic equation which has the solution 



1 / fd? ~ m x Q x\ — Xq ■ xi d\ ( ' m 2 \ 

^ = 2(r +4(1+7) ^ r Xi (15) 

The solution with the upper (lower) sign is denoted by ip\ (^2)- For astrometry the 
solution <f 1 can be considered to be the more relevant solution, because y>2 represents 

(m 2 \ 
~d/ 2 ) 

are smaller than or equal to — Equation (TT51) represents the generalized lens 

4 d' 

equation. This equation is applicable not only for any configurations where the post- 
Newtonian approach, the post-post-Newtonian approach, or the classical lens equation 
is valid, but also for all those extreme configurations given by the points (1) - (3) in the 
introductory section. Especially, it allows an analytical investigation of light-deflection 
in binary systems [20J. In the following Section we will show that the formula ( TToT) 
represents a link between standard post-Newtonian approach and classical lens equation. 

5. Discussion of generalized lens equation 

5.1. Comparison with standard post-Newtonian and post-post-Newtonian approach 

In this Section we compare the generalized lens equation ( TT5l) with the standard post- 
Newtonian and post-post-Newtonian approach of light-deflection. A series expansion of 
the solution if\ in Eq. ( TT5T) for d^> m yields 



Vi = <^ P N + ^ PP N + O (m 3 ) + O f ^ j > 



(16) 



with 



<^ P N =(1 + 7) -7 p - 4 T' 

a Hxi d 

2 m 2 {xqXi - x ■ x x ) 2 m 2 x x 

^= - (1 + 7) ~¥ Wdx^ - 16 ^T' (18) 

Expression (TTTT) is called standard post-Newtonian solution, cf. Eq. (24) in [UJ. The 
expression (TTgj) is just the 'enhanced' post-post-Newtonian term, cf. Eqs. (3) and (4) 
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in [rj]. The 'enhanced term' can be arbitrarily large for small d and large x\. That 
is the reason why the standard post-Newtonian and post-post-Newtonian solution is 
not applicable for extreme configurations like binary stars. The term O (m 3 ) will be 
discussed below, see Eq. ([27]) ; here it is only essential to realize that this term can be 



larger than the neglected terms O — 
6 6 ' d' 



m 



5.2. Comparison of generalized lens equation and classical lens equation 

The standard weak-field lens equation is usually called classical lens equation and given, 
for instance, in Eq. (67) in [9] or Eq. (24) in |JJ. Let us briefly reconsider the classical 
lens equation. According to the scheme in Figure [2] we obtain the following geometrical 
relations 

<P + ip = 5, (19) 
a tan <p = b tan if) . (20) 

Here, the angles are if> = 8(fi, k) and 5 = 5(n, fjt), where /x = is the unit tangent 

vector at the position of the source in the direction of the propagation of the light 
signal. If the source and observer are infinitely far from the massive body, then the 
total light-deflection angle 5 = 5 (n, /x) in Schwarzschild metric reads [3] 



i = 2(l +7 )- + (- p ) 



(21) 



j can be estimated 



to be smaller than or equal to — ^ , see [3| . In classical lens approach, the 



15 Ti m 2 

approximation d' ~ d + a tany? is used, see Figure [2 Inserting this relation into 
f l2T|) . by means of geometrical relations f fT9|) and (|2T)]). and using tan ip = ip + 0(p 3 ) and 
temip = if) + 0((p 3 ), we obtain the quadratic equation 

n d , s m b , . 

^ + -^-21 + T -- - = 0. (22) 
a a a + b 

The solution of Eq. ( 122]) is the classical lens equation: 




a 2 




rfr = ;h/^ + 8(l + 7)-— tT-I, (23) 



which is valid in case of a, b ^> d; the solution with the upper (lower) sign is denoted by 

„class I „class\ 

It should be noticed that in (125]) not only the light deflection angle (p is assumed to 
be small, but also the source and observer are assumed to be far from the massive body, 
i.e. 5(x , xi) ~ 7r; note that due to that fact equation ( 123]) agrees with the classical lens 
equation (67) in (9]. Therefore, the classical lens equation is not applicable for extreme 
configurations like binary systems or light deflection of moons at their giant planets of 
solar system. 
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It can easily be shown that the classical lens equation fl23l) follows straightforward 
from the generalized lens equation (fT5l) . That means, if we rewrite (fT5l) in terms of 
a = k ■ X\ and b = —k ■ xq and perform a corresponding series expansion of generalized 
lens equation ( TT5|) then we just obtain the classical lens equation (|23|) as the leading 
term in this series. 

Furthermore, in the limit d — > 0, known as Einstein ring solution, the generalized 
lens equation ( TT5l) and the classical lens equation ( 1231) yield the same result: 



lim = lim < 2 SS = 4/2 (1 + 7 ) - ■ (24) 

Finally, we note that in the extreme configuration b = (in this limit <f2 does not exist) 
we obtain from (I15p the result 



lim ip 



d 2 , , m da d \ \, . m 



while the classical lens equation yields simply Lpf ass = 0. Obviously, in the limit a — > 
the expression (|25|) yields zero as it has to be because in this limit the distance between 
source and observer vanishes, that means no light deflection. 

5.3. Comparison with exact solution 



The accuracy of (1151) and the stated estimate that the neglected terms are smaller than 
or equal to — -j-^ has also been confirmed by a comparison with the exact numerical 
solution of ©. 

For that, we have solved the geodesic equation ([2]) in Schwarzschild metric by 
numerical integrator ODEX [10] for several extreme astrometrical configurations. Using 
forth and back integration a numerical accuracy of at least 10 -24 in the components of 
position and velocity of the photon is guaranteed. Thus, the numerical integration can 
be considered as an exact solution of geodesic equation which is denoted by ^ num . This 
numerical approach has been described in some detail in [TT] . In all considered extreme 
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Figure 3. Comparison of solution <p\ of generalized lens equation (ITSt with 
exact numerical solution (p n um for the case of a grazing ray at Sun (A) (x% = 
(-la.u.,O,O),m = 1476.6 m, d' = 696.0 x 10 6 m) and Jupiter (B) {x 1 = 
(-6.0a.u.,0,0),m^ = 1.40987m, d' = 71.492 x 10 6 m), where a.u. = 1.496 x 10 11 m 
denotes astronomical unit. 



configurations the validity of ( TT5]) and the given estimate of neglected terms have been 
confirmed. As example, in Figure [3] we present the results for light-deflection of a 
grazing ray at Sun and Jupiter. These examples elucidate the fact that the accuracy of 
generalized lens equation ( Il~5l) is much beyond microarcsecond level of accuracy in case of 

light-deflection at giant planets. The reason for this fact is that - <C //as for giant 

planets of the solar system; only in the vicinity of the Sun we have - — g ~ 11 £*as. 

The accuracy shown in Figure [3] (B) is considerably better than the post-post- 
Newtonian solution investigated in detail in [TT] US] , cf. Figure |3] (B) with FIG. 2 in 
[18J. In order to understand the numerical difference between Figure [3] (B) and FIG. 2 
in [18], we perform a further series expansion of Eq. ( I15p up to terms of order m 4 , that 
means 



(26) 



Vl = ¥>pN + ^ppN + ^PPPN + O (m 4 ) + O ^^72 j 
where the 'enhanced' terms beyond post-post-Newtonian terms are: 

^ pppN = 2(l + 7 ) <128^^. (27) 

The given estimation in ( 1271) shows, that for large x\ this term can be considerably 

(vr? \ 
J ■ Moreover, the numerical difference 

between Figure |3] and FIG. 2 in [18] is just given by the term (j2~ 



6. Summary 

Modern astrometry has achieved a microarcsecond level of accuracy, e.g. astrometric 
missions SIM (NASA) or GAIA (ESA). A time-efficient data reduction implies the need 
of approximative and highly precise solutions for the light deflection on this level of 
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accuracy. In our investigation we have suggested a generalized lens equation ([15]) for 
weak gravitational fields of Schwarzschild metric and valid for finite distances of source 
and observer from the light deflecting body. The derivation is based on the solution 
of geodesic equation (fTTj) in post-Newtonian metric and Chandrasekhar's coordinate 
independent impact parameter d' fl6]) and its relation to the light-deflection angle ip 
given in f JT4|) . The neglected terms in ( T15]) can be estimated to be smaller than or 

equal to — ^— The accuracy of generalized lens equation (ITS]) is considerably better 
than the standard post-Newtonian and post-post-Newtonian approach, which has been 
investigated in some detail in [TTJ [18] and the reason for this fact has been pointed out. 

The generalized lens equation (I15p satisfies three conditions formulated in the 
introductory Section. Moreover, we have shown that in the appropriate limits we obtain 
the post-Newtonian terms, 'enhanced' post-post-Newtonian terms and the classical lens 
equation. Thus, the generalized lens equation ( fT5l) provides also a link between these 
essential approaches to determine the light-deflection. Numerical investigations have 
confirmed the analytical results obtained. 

The generalized lens equation ffl5|) allows an analytical understanding and 
investigation of light-deflection in extreme astrometric configurations. Especially, the 
determination of light-deflection in binary systems using of generalized lens equation 
( ]1~5"]) has been investigated in [20] . 
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